Fast DNA translocation through a solid-state nanopore 
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We report translocation experiments on double-strand DNA through a silicon oxide nanopore. 
Samples containing DNA fragments with seven different lengths between 2000 to 96000 basepairs 
have been electrophoretically driven through a 10 nm pore. We find a power-law scaling of the 
translocation time versus length, with an exponent of 1.26 ± 0.07. This behavior is qualitatively 
different from the linear behavior observed in similar experiments performed with protein pores. 
We address the observed nonlinear scaling in a theoretical model that describes experiments where 
hydrodynamic drag on the section of the polymer outside the pore is the dominant force counter- 
acting the driving. We show that this is the case in our experiments and derive a power-law scaling 
with an exponent of 1.18, in excellent agreement with our data. 

PACS numbers: 8714Gg,8715Tt 



Translocation of biopolymers such as polypeptides, 
DNA, and RNA is an important process in biology. Tran- 
scribed mRNA molecules for example are transported out 
of the nucleus through a nuclear pore complex. Viral 
injection of DNA into a host cell is another example. 
Translocation of DNA and RNA can be studied in vitro, 
as demonstrated by Kasianowicz et al. jl| using an a- 
hemolysin pore in a lipid membrane. By measuring the 
ionic current through a voltage-biased nanopore, one can 
detect individual single-strand molecules that are pulled 
through the pore by the electric field. Li et al. 0, El 
showed that solid-state nanopores can also be used for 
such experiments. We report a set of experiments with 
silicon oxide nanopores on double-strand DNA with vari- 
ous lengths. Surprisingly, we find a nonlinear scaling be- 
tween the translocation time r and the polymer length 
Lq, in contrast to the linear behavior observed for all ex- 
periments on a-hemolysin 0, • In our experiments we 
find a clear power-law relation r ~ L^ 26 , for DNA frag- 
ments from 2000 to 96000 basepairs (bp). While a com- 
plete model for translocation should in principle include 
hydrodynamic, steric, electrostatic and entropic effects, 
we argue that in our experiments the dominant contri- 
butions to the force balance come from hydrodynamics 
and driving, and propose a simple model that accurately 
reproduces the observed scaling. 

Experimental results. Figure \B(a) shows the experi- 
mental layout for translocation studies. At the heart of 
the setup is a solid-state nanopore device, fabricated by 
shrinking a 20-50 nm pore in silicon oxide in a trans- 
mission electron microscope to a final size of 10 nm 
[|| . The nanopore is situated in an insulating membrane 
which separates two macroscopic reservoirs filled with an 
aqueous buffer solution. When a voltage bias is applied 
over the membrane in the presence of DNA molecules 
in the negative compartment, the DNA is electrophoret- 
ically drawn through the pore due to its negative charge. 
The detection technique is simple and elegant: A poly- 
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FIG. 1: (a) Cross-sectional view of the experimental setup 
(not to scale). A negatively-charged DNA molecule is elec- 
trophoretically driven through a 10 nm aperture in a silicon 
device, which is located between two reservoirs kept at a po- 
tential difference. Both reservoirs are filled with an aque- 
ous buffer solution (1 M KC1, 10 mM Tris-HCl pH 8.0, ImM 
EDTA). (b) Measured ionic current versus time, after the ad- 
dition of 4 /im (11.5 kbp) DNA to the top reservoir. An 
individual event is shown with increased time resolution, (c) 
Dwell time vs. polymer length. The line shows the result 
of a power-law fit to the data, with a best-fit exponent of 
a = 1.26 ±0.07. 



mer traversing the pore lowers the amount of conduct- 
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ing solution present inside the pore and thereby reduces 
the ionic conductivity between the reservoirs. Passing 
molecules are thus detected as short dips in the ionic 
current, which is induced by the externally applied volt- 
age (see Fig. [QffcJ). Analogous to Li et al. we find 
that the molecules with a diameter of about 2 nm can 
pass the 10 nm pore either in a linear or in a folded 
fashion. Using an event sorting algorithm discussed else- 
where , we analyze only the linear unfolded transloca- 
tion events in this work. Fig. GJf6J shows an example of 
such a linear translocation event, detected with 11.5 kbp 
linear DNA. The width of the dip is interpreted as the 
duration of the translocation. We performed three exper- 
iments, all at room temperature: One on linear 11.5 kbp 
DNA, one on linear 48 kbp A-DNA (here we detected 
both individual molecules and dimcrs of two molecules 
bound together with their complementary sticky ends), 
and one on a mixture that contains 27491 bp, 9416 bp, 
6557 bp, 2322 bp, and 2027 bp fragments. In the last 
experiment, the length difference between the 2322 and 
2027 bp fragments could not be resolved experimentally. 
The durations of individual linear events were collected 
in a histogram for each experiment. Figure [T]fc ) shows 
the dwell time (determined as the peak position in these 
histograms) versus polymer length for all 7 DNA frag- 
ments. The full width at half maximum of the peaks are 
taken as error bars. We find a clear power-law scaling 
of the dwell-time r with the length Lq, t ~ Lq , with a 
value of 1.26 ± 0.07 for the exponent a. These experi- 
ments are discussed in more detail in |6J . It is interesting 
to note that our results are in good agreement with those 
reported by Li et al. 3]. They have studied the translo- 
cation of 3000 bp and 10000 bp DNA through a 10 nm 
silicon nitride nanopore, and find translocation times of 
100 /is and 400 /is respectively. 

The translocation process consists of two separate 
stages. First, there is the capture stage. A DNA molecule 
initially in solution in the negative reservoir has to come 
close enough to the pore to experience the electrostatic 
force and get pulled in. We assume that the reservoirs 
are good ionic conductors, and the driving force is only 
felt in the direct vicinity of the pore. Capture is thus a 
stochastic process, since the pore has to be reached by dif- 
fusion. In this work, we focus on the second stage, where 
the DNA passes the pore until it has reached the other 
side. We assume that one end of the DNA has entered 
the pore and calculate the time required for complete 
translocation. 

Slow vs. fast translocations. We now address the de- 
pendence of this duration on the length of the polymer. 
To this end, we consider a linear polymer consisting of N 
monomers, each of which has a Kuhn length b. This poly- 
mer is partially threaded through a narrow pore. Time 
t = sets the moment of initial capture. We will let L(t) 
denote the contour length of the untranslocated part of 
the polymer, so that L(0) = Nb = Lq. The dwell time 



r is therefore determined by L(t) — 0. A second time 
scale in the problem is the characteristic relaxation time 
scale of the translocating polymer. This Zimm time Q, 
given approximately by 

riR s n 

t z *0A-^, (1) 

can be considered an upper bound on the time it takes 
the polymer to relax to an entropically and sterically fa- 
vored configuration. In this expression, rj is the solvent 
viscosity and R g is the radius of gyration of the polymer. 
This is the radius of the typical blob-like configuration 
that the polymer will assume in a good solvent, and it 
scales with the polymer length as 

Rg ~ Lq, (2) 

which defines the swelling exponent v. It depends on the 
dimensionality of the system, and theoretically a value 
of 0.588 is found for self-avoiding polymers in a good 
solvent ||. Smith et al. @ have measured the diffu- 
sion constant D for stained DNA molecules with lengths 
ranging from 4.3 kbp to 300 kbp. They report a scal- 
ing with length L as D ~ L~ v with v = 0.611 ± 0.016, 
and conclude that Flory scaling is appropriate for DNA 
molecules longer than about 4 kbp. At room tempera- 
ture, the measured velocity is about 0.8 /us per base or 
slower A 100-base, single-stranded DNA fragment 

therefore takes around 80 fis to fully translocate. When 
we compare this to the Zimm time for the same poly- 
mer fragment, about 0.2 /xs, we see that relaxation is 
much quicker than the translocation. We will call such 
events, for which t > tz, slow translocations. Luben- 
sky and Nelson have argued that for single-stranded 
DNA and RNA through a-hemolysin, the criterion for 
slow translocation is indeed satisfied for polymer lengths 
up to hundreds of nucleotides. They show that the Zimm 
time for a polynucleotide of roughly 300 bases is compa- 
rable to the translocation time per nucleotide. 

The criterion for slow translocation is evidently not 
met in our experiments on solid-state nanopores. A full 
A-phage genome (48.5 kbp, or 16.5 /im of double-stranded 
DNA) is found to take only around 2 ms to traverse a 10 
nm SiC>2 pore. The Zimm time for this molecule, in com- 
parison, is about 700 ms, clearly much longer than the 
translocation time. Even the translocation of the short- 
est molecules studied in these assays (2000 bp) can be 
considered fast, with a translocation time of about 50 /xs 
and a Zimm time of just over 2.3 ms. We therefore refer 
to this second regime, where r <C tz, as fast translo- 
cations. We should point out that an important reason 
for the fastness of our system is the fact that we use 
double-stranded DNA, which has a much larger persis- 
tence length than single-stranded DNA and consequently 
has a long relaxation time. 

Model. Let us estimate the magnitudes of the possi- 
bly relevant forces, following Lubensky and Nelson [ll| . 
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FIG. 2: At time t the DNA on the left side of the pore has 
a radius of gyration of R g , indicated by the dashed line. The 
balance between the two dominating forces determines the 
dynamics: A driving force that locally pulls the DNA through 
the pore and a viscous drag force that acts on the entire DNA 
blob. 

First, consider the driving. As stated, a potential dif- 
ference across the pore exerts a highly localized force 
on the negatively charged DNA molecule. We assume 
the potential drop to occur entirely inside the pore, and 
therefore only the part of the polymer inside experiences 
the driving force. This force can then be estimated as 
-^driving = 2eV/a, where e is the elementary charge, V 
is the potential difference and a = 0.34 nm is the spac- 
ing between nucleotides. A bias voltage of 120 mV, as 
is typically used in experiments, thus produces a force 
of around 110 pN. This value is an upper bound of the 
actual force, as screening effects may greatly reduce the 
effective charge on the DNA, and thereby the driving 
force. Simulations of Manning condensation on double- 
stranded DNA yield charge reduction values between 
53% and 85% [HE1- Barring complete screening how- 
ever, we consider our DNA translocations to be strongly 
driven and this justifies ignoring diffusive contributions. 

In the absence of specific DNA-pore interactions, the 
viscous drag per unit length in the pore can be estimated 
as 2irr]rv/(R — r), where R is the pore radius, r is the 
polymer cross-sectional area, rj the solvent viscosity and 
v is the linear velocity of the polymer inside the pore. 
Substituting typical values (77 = 1.10 -3 Pa-s, r = 2 nm, 
v = 10 mm/s, R — 10 nm and a pore depth £ porc of 
20 nm) we can estimate this drag force to be around 
0.3 pN, decidedly smaller than the driving force. We feel 
this constitutes an essential difference between solid-state 
pores and protein pores: In sufficiently shallow solid state 
pores the effect of friction inside the pore is negligible. 

Finally we estimate the hydrodynamic drag on the un- 
translocated part of the polymer, outside the pore. To 
this end, we approximate the untranslocated part as a 
sphere of radius R g (see Fig. [2J - As the polymer threads 
through the pore, the center of mass of this sphere moves 
towards the pore at a velocity dR g /dt. Assuming Zimm 
dynamics (and thus that the solvent inside the coil moves 



with the polymer), the coil experiences a Stokes drag 
force of Qit-qRv — 6irrjR g AR g / 'At, which for typical pa- 
rameters yields a drag force of about 24 pN. This as- 
sumption is justified by experiments by Smith et al. 0], 
who found clear evidence for Zimm dynamics for DNA 
longer than 4.3 kbp. Clearly, in this case the hydrody- 
namic friction on the part of the polymer outside the pore 
is the dominant force counteracting the driving force. We 
therefore choose to model fast translocation dynamics as 
determined only by the cumulative effect of driving at 
the pore and hydrodynamic friction outside. 

FigureEldepicts the simplified system we consider. The 
part of the polymer inside the pore experiences a driving 
force to the right, while the length of polymer before the 
pore is coiled up. The pore is sufficiently small to allow 
only linear (i.e. unfolded) passage of a single molecule at 
a time. 

As the polymer is pulled through the pore the blob 
before the entrance shrinks in size, and thus its center of 
mass moves towards the pore with a velocity 

^blob = R g ~ L V ~ X L , (3) 

where the dot denotes a time derivative. Motivated by 
our consideration of the relative magnitudes of the coun- 
teracting forces, we propose the principal effect of hy- 
drodynamics is to resist motion with a Stokes drag on 
the DNA coil that is proportional to its radius times the 
velocity, 

i"drag ~ -RVblob ~ RgRg 

~ L 2v ~ x L . (4) 

Force balance must be met at all times, and since there 
are only two major forces the driving force should bal- 
ance the hydrodynamic friction: i*d rag = —F, driving- As 
the driving force is constant during the whole process, 
the same holds for -Fd ra g. Thus we can extract the lin- 
ear velocity vy ln = —L of the DNA inside the pore from 
Eq. |gj): 

v lin = -£ ~ -L 1 ^, (5) 
which allows us to obtain the dwell time r by integration, 

r= f T At= [ v^(L)AL~Ll». 
Jo JL 

On the basis of this model we thus predict a powerlaw 
relation between the dwell time r and the contour length 
Lq. Taking the theoretical value of 0.588 for the Flory 
exponent v we find an exponent of a = 1v = 1.18 for this 
model. If we take the experimentally obtained value for 
v of 0.61 0, we find a — 1.22 - in excellent agreement 
with our experiments, where we find power law scaling 
with an exponent of 1.26 ± 0.07. 
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Scaling regimes for translocation. General considera- 
tions along the lines of the argument presented in the pre- 
ceding sections can be used to qualitatively understand 
the various regimes of polymer translocation. Firstly, 
it is important to determine the dominant contribu- 
tion to the friction. In most cases, it suffices to com- 
pare the pore friction F porc = £ e ffUu n (with £ e ff equal 
to 2irr]lp OTC r/(R — r) in the absence of specific interac- 
tions) to the Stokes drag on the coil dirriRgRg. If the 
pore friction dominates, force balance with respect to 
the constant driving force implies that the translocation 
time scales linearly with the polymer's length r ~ L . 
A possible reason for a large pore friction could be the 
presence of specific interactions, but because of the ge- 
ometric factor in the effective friction constant £ e ff the 
shape of the pore could also lead to pore friction domi- 
nated translocation. Such linear dependence of r on the 
length for single-stranded DNA ranging from 12 to 400 
bases has been reported experimentally by Kasianowicz 
and Meller Q . For the a-hemolysin pore they used it 
is indeed speculated that significant specific interactions 
with the passing DNA occur. 

When the Stokes drag dominates one can derive, with- 
out any assumptions on the polymer statistics, that 
t ~ R 2 g . Depending on the length of the polymer dif- 
ferent regimes are thus obtained: when the polymer is 
short compared to its persistence length R g ~ Lq, and 
we find that r ~ Lq. For polymers of intermediate length 
the radius of gyration follows the scaling for a Gaus- 

1 12 

sian chain, R g ~ L , and consequently the translo- 
cation time is predicted once again to scale linearly with 
length (note, however, that this is a qualitatively different 
regime than the pore- friction dominated regime identified 
before). For long polymers (such as those considered in 
the preceding sections) we have shown that r ~ Lq v . 

So far we have presented scaling arguments assuming 
Zimm (non-free draining) dynamics. Kantor and Kardar 
have identified and numerically confirmed yet another 
regime [lH ] where a = v + 1. This behavior can be un- 
derstood assuming Rouse dynamics (stationary solvent). 
In this case, the hydrodynamic drag would be given by 
-Fdrag ~ Lvuob , and one does indeed recover their scaling 
law t ~ Lq +v '. We speculate that this regime might be 
observable for semidilute solutions close to c*. 

In all cases considered in this section, one can indepen- 
dently determine whether or not the polymer is frozen in 
its configuration by comparing the Zimm- and translo- 
cation timescales. We do not expect this to affect the 
scaling behavior in the various regimes, but it will affect 
the prefactors. 

Concluding remarks. We have obtained a simple and 
elegant model description that appears to describe our 
data well. There are however several effects we neglect 
but which could have an additional influence on the pro- 
cess that we consider. For instance, we ignore any fric- 



tion experienced by the DNA that has already passed 
the pore. We also expect that an electro-osmotic flow 
is generated inside the pore. This effect is caused by an 
electrophoretic force on the ions screening the charge on 
the surface of our pore. As silicon oxide is known to be 
negatively charged in water, there is a surplus of posi- 
tive ions near the surface. These positive ions generate a 
flow of water inside the pore, slowing down the DNA that 
moves in the other direction. While we have not explored 
the consequences of these possibilities the observed agree- 
ment between theory and experiment suggests that at 
least for the fast polymer translocations considered here, 
hydrodynamic drag does indeed dominate the dynamics. 
Identification and understanding of the dominant effects 
in polymer translocation through nanopores is relevant 
not only for biological processes, but also for potential an- 
alytical techniques based on nanopores. Rapid oligonu- 
cleotide discrimination on the single-molecule level has 
been demonstrated with a hemolysin and more re- 
cently solid-state nanopores were used to study folding 
effects in double-stranded DNA molecules |3J, 111 • Future 
applications of this technique may include DNA size de- 
termination, haplotyping and sequencing. 
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